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ABSTRACT 

Galaxies evolve continuously under the influence of self-gravity, rotation, accretion, mergers and 
feedback. The currently favored cold dark matter cosmological framework, suggests a hierarchical 
process of galaxy formation, wherein the present properties of galaxies are decided by their individual 
histories of being assembled from smaller pieces. However, recent studies have uncovered surprising 
correlations among the properties of galaxies, to the extent of forming a one-parameter set lying on 
a single fundamental line. It has been argued in the literature that such simplicity is hard to explain 
within the paradigm of hierarchical galaxy mergers. One of the puzzling results, is the simple linear 
correlation between the neutral hydrogen mass and the surface area, implying that widely different 
galaxies share very similar neutral hydrogen surface densities. In this work we show that sclf-rcgulated 
star formation, driven by the competition between gravitational instabilities and mechanical feedback 
from supernovae, can explain the nearly constant neutral hydrogen surface density across galaxies. We 
therefore recover the simple scaling relation observed between the neutral hydrogen mass and surface 
area. This result furthers our understanding of the surprising simplicity in the observed properties of 
diverse galaxies. 



Subject headings: galaxies: kinematics and dynamics 
galaxies: stellar content 



galaxies: ISM — supernova remnants 



1. INTRODUCTION 

In a hierarc hical process of galaxy formation 
(|Cole et al.ll2000[) . within the currently favored A-CDM 
cosmological framework, the present properties of galax- 
ies are decided by their indi vidual h istories of being as- 
sembled from smaller halos (iLaccy fc Cole 1993)- More- 
over, these galaxies evolve continuously under the in- 
fluence of self-gravity, r o tation , accretion and feedback. 
However, iDisnev et al.l ()2008D have recently reported 
surprising correlations among the properties of galax- 
ies, to the extent of forming a one-parameter set ly- 
ing o n a single fund amental line (j Garcia- Appadoo et al.l 
2009( 1. Ivan den Bergh (,2008[ ) has argued that such sim- 
plicity is hard to explain within the paradigm of hier- 
archical galaxy mergers. Qpe of the puzzling results 
(iGiovanelli fc Havnesi 119831: iHavnes fc Giovanellil 1 1984 : 



Verheiien & Sancisi"2001': ' Rosenberg fc Schneideil 12003 
Garcia- Appadoo ct al. 2009), is that the neutral hydro- 
gen mass across these galaxies scales almost linearly with 
surface area, implying that widely different galaxies share 
very similar neutral hydrogen surface densities. This 
nearly constant HI surface density has been pointed out 
in literature to be an intriguing puzzle, demanding an 
explanation. 

In this work wc argue that correlation between the neu- 
tral hydrogen mass and surface area may be preserved, 
despite the complex merger histories of the galaxies, by 
self-regulated star formation, driven by the competition 
between gravitational instabilities in the rotating disk 
and mechanical feedback from supernovae. When merg- 
ers drive a galaxy away from the fundamental line, self 
regulation of the porosity of the ISM and gravitati onal 
instability of the star forming disk, as proposed by iSilkl 
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can bring the neutral hydrogen surface density 
back to the value which is predicted in our simple model. 
This can explain the regulation of the neutral hydrogen 
surface density in galaxies, explaining part of the surpris- 
ing simplicity in the observed properties of galaxies. 

2. SURFACE DENSITY OF GAS IN GALAXIES 

Evidence of a nearly universal neutral hydrogen (HI) 
surface density in galaxies has been accumulating over 
the past three decades. The initial hints came from 
single-dish 21 cm radio obs ervations of nearby galaxies. 
IGiovanelli fc Havned (jl983f ) reported neutral hydrogen 
(HI) observations of 24 galaxies in the Virgo cluster, us- 
ing the Arecibo telescope. The HI sizes and masses were 
found to be correlated in the same manner, irrespective 
of whether they were HI rich or HI poor. A similar cor- 
relation was soon found between the optical sizes and 
HI m asses of 288 isolated galaxies (jHavnes &: Giovanellil 
IT981 . It was also found that the optical diame- 
ters of spiral disk are better cor related with the HI 
mass than the morphologica l tvpe (IHavnes fc G iovancflil 
I1984D . A deeper survey bv I Verheiien fc Sancisil (|200!l ) 
has revealed similar correlations for galaxies in the 
Ursa Major Cluster. The A recibo Dual-Beam Survey 
(|R,osenberg fc SchneideH (200(1 (ADBS) has found HI in 
265 galaxies in a "blind" survey of ^ 430deg^ of sky. 
While most of the ADBS galaxies were unresolved at the 
resolution of the Arecibo, the Very Large Array (VLA) 
was used for interferometric mapping of 84 galaxies and 
determine accurate sizes of 50 of them. A comparison 
of HI masses and HI sizes of these along with 53 galax- 
ies with high resolution maps from literature revealed 
that they were consistent with a nearly constant av- 
erage HI surfac e density of the order ^ 1O''M0 kpc~^ 
(jRosenberg fc Schneider. .20031) . The HI masses of indi- 
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vidual ADBS galaxies, spanning 3 orders of magnitude, 
deviate by only ^ 0.13 dex (Icr) from those expected 
from a constant HI surface density. This puts the evi- 
dence, for a regulated HI surface density across galaxies, 
on a firm observational basis. 

A recent study of HI selected galaxies (free from op- 
tical selection effects) found using the Parkes radio tele- 
scope and identified with SDSS sources, has shown that 6 
observed parameters, namely the dynamical mass (Md), 
HI mass (Mhi), luminosity, color, and two radii contain- 
ing 50% and 90% of the observed luminosi ty, have 5 in- 
depen dent correlations among themselves (jPisnev et al.l 
l2008f ). This implies that the galaxies form a single pa- 
rameter family and are not removed significantly from 
their fundamental line by the diverse merger histo- 
ries that they would have had in the process of hi- 
erarchical galaxy formation. Some of the correlations 
have already been widely know n and discussed in o ther 
forms, such as the correlation (jGavazzi et al.lll996D be- 
tween luminosity and dynam i cal m ass. iDisnev et al.l 
(|2008[ ): iGarcia-Appadoo et al.l (|2009D report a tight lin- 
ear correlation in the already established relation be- 
tween the HI mass and the surface area. The nearly 
constant HI surface density is argued in literature to 
be an intriguing puzzle whic h demands an explanation 
(|Garcia-Appadoo et al.ll20d9l ). We discuss below, a phys- 
ically motivated explanation for this important observa- 
tion. 

3. GRAVITATIONAL INSTABILITY IN DISK GALAXIES 

Below a surface density threshold, azimuthally inte- 
grated star form ation in giant HII regions across a galaxy 
ceases (jKennicut t 1989.1. The existence of this threshold 
is traditionally stated in the for m of the Toom re param- 
eter for gravitational instability ()Spitzeiill968() . Below a 
Toomre parameter Q < 1, rotation support and gas pres- 
sure cannot stabilize a thin self-gravitating disk against 
gravitational instabilitv .Safronov. (,1960.) : ,ToomrG(,1964 . 
Here Q is defined as 



ity as 



gas 5 



(1) 



where the critical gas surface density is given by fXcr = 
flag/nG, in terms of the angular ve locity fl and v eloc- 
ity dispersion ag in the gas disk. iSchavd ()2004[ ) has 
demonstrated that results from detailed considerations 
including not only self-gravity, but metals, dust and UV 
radiation, coincide with this empirically derived surface 
density threshold for star formation. 

This condition depends on the local angular velocity, 
which we may define the as = V x V. Using cylin- 
drical polar coordinates we can now write down the curl 
operator as 



V X A = - 
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Assuming, rotation and translation symmetry along (f) 
and z respectively, if the bulk motion of the gas is only 
in the (f) direction (as would be the case if the motion is 
purely Keplarian) we may write down the angular vcloc- 
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where the radial coordinate p has been replaced by r, to 
avoid confusion with density. In nearly rigid rotors such 
as dwarf galaxies, the velocity is given by V(r) = i^pr. 
In su ch a case, | V x V| is given by (lArfken fc Weberi 
I2OOI problem 2.4.7). For flat {V{r) = Vq) rotation 
curves |V x V| ~ Vq/t. Hence, for all these cases, at 
boundary of the starforming disk, the local value of 
the angular velo city i s com parable to the global value. 
Hence, following |Sil3 (Il997l) we shall use the global an- 
gular velocity of the disk in the rest of this work. 
The global angular velocity can now be expressed as 



GMd 
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in terms of the enclosed dynamical mass within the gas 
disk. Replacing Md = AnpdR^/S, w here Pd is the nearly 
universal dynamical mass density (iDisnev et al.l 120081: 
IGarcia-Appadoo et al.ll200"9() . we have flo = y/AirGpd/i- 
Substituting, we get 
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Now, we note that for the star forming disk to be long 
lived, it may only be marginally unstable. Q is seen to 
be 1 throughout t he star forming regions in well ob- 
served disk galaxies (|Kennicuttlll989f ) . This implies that 
throughout the star forming disk the mean gas surface 
density is of order pgas ~ p-cr- Hence, the mean surface 
density of gas is controlled primarily by the dynamical 
mass density pd and the gas velocity dispe rsion an- 

One of the tight correlations see n by IDisnev et al.l 
(|2008f)HGarcia-Appadoo et al.l (I2009D is between the dy- 
namical mass and cube of the optical radius, implying a 
roughly constant average dynamical mass density (within 
the radii of their star forming disks) for galaxies. It is ex- 
pected that Dark Matter contributes most of the gravita- 
tional mass of the galaxies, hence a physical explanation 
of this observation would require an understanding of the 
nature of Dark Matter. Whether or not hierarchical halo 
build-up within a A-CDM cosmology c an explain this 
observation is still under investigation. iLoeb fc Weined 
(i20ld) have suggested that cold dark matter particles 
interacting through a Yukawa potential could make ha- 
los beyond a certain critical density evaporate over an 
Hubble time. This may set a characteristic scale to the 
peak density of dark matter halos. Note however that 
the dark matter halos of galaxies are much larger than 
the sizes of their star forming disks. As a result, most of 
the dark matter may lie further out. iKentI (jl987| ) points 
out that the relative contributions of the dark matter 
halo and stellar contents to the dynamical mass of a 
galaxy vary signiflcantly with its luminosity and mor- 
phological type. In the rest of this work, we shall use 
the simple Md (X- relation observ ed by IDisnev et al.l 
(|2008| ): IGarcia-Appadoo et all (|2009[ ). showing a shared 
Pd ~ 10^ MQkpc~^ across galaxies. This implies, that the 
mean surface density is controlled essentially by the gas 
velocity dispersion ag which is driven by energy input 
from supernova explosions in the disk. 
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4. MECHANICAL FEEDBACK FROM SUPERNOVAE 

The distribution of gas in our galaxy has been vari- 
ously describe d in the past as bei ng similar to that of 
Swiss Che ese (ICox fc SmithI [l97l . as a Cosmic Bub- 
ble Bath (iBrand fc Zealevlll975D or as the Violent In- 
terstellar Medium (jMcCrav fc Snowlll979D . The basic 
idea is that supernova remnants are full of coronal gas 
which can persist for long enough, as it radiates very 
inefficiently by bremsstrahlung, so that a modest su- 
pernova rate may produce an interconnec ted morphol- 
ogy of hot gas. Shells and supershells (|Heilej |1979() 
have been fo und in the neutral hydro gen distribution 
of galaxies. iMcCrav fc KafatosI (|1987f ) suggested that 
stellar winds and repeated supernovae from an OB as- 
sociation may create cavities of coronal gas in the inter- 
stellar medium lead i ng to the formation of supershells. 
iChakraborti fc Ravi ()201lD have demonstrated that this 
mechanism can explain the dynamics of a HI supcrshell 
found in MlOl driven by mechanical energy input from 
supernovae in a giant young stellar association. 

The fraction of volume filled by supernova-driven hot 
gas is referred to as the porosity (P) of the interstel- 
lar matter (ISM). The porosity is driven by the supply 
of hot gas from supernova remnants (SNRs) and is re- 
lated to th e 4-volume of a SNR in the cooling phase 
(jCioffi et al . 1988) (i'sn) and the supernova rate per unit 
volume (tsn) as 

P = usN X rsN- (6) 

The supernova rate is related to the recent star formation 
rate per unit volume (p*) as 



rsN P*/'msN, 



(7) 



where msN is the total mass of star formation required 
on an aver age to produce each core collapse supernova. 
iSilkl ()1997[ ) has argued that if P is too high, it would 
suppress the efficiency of star formation. Since the su- 
pernova rate follows the recent star formation rate, the 
situation called "blowout" would throttle the supply of 
hot gas and bring down the value of P. On the other 
hand if P was too low, it would allow the cold gas phase 
to dominate and form new stars more efficiently, some 
of which would soon explode as supernovae and drive up 
the value of P. Hence, there would be a self regulation 
process that woul d con t rol P. 

In this manner, ISilS (|1997[ ) shows that supernova ex- 
plosions supply the energy input necessary to maintain 
the velocity dispersion in the gas phase at 



(Ta = 6.90P 



-0.58^0.1 ^0^2 ^0.008]^jj^g- 



(8) 



This expression depends only weakly on the gas density 
Ug, energy input from individual supernovae £'51 x 10^^ 
ergs, and the metallicity Q all of which have typical val- 
ues of order unity. The dependence on these parameters 
is dropped from our subsequent equations. Assuming 
self-regulated star formation [P ^ 0.5) the predic ted 
gas velocity dispersion has been shown (jSilkl 119971 ) to 
be an ^ 11 km s~^, close to the value of observed by 
iStark fc Brandl ()1989[ ) for the three-dimensional peculiar 
velocity dispersion of interstellar molecular clouds within 

3 kpc of the S un. 

iJoung et al.l ()2009[) have studied the effect of the su- 
pernova rate on the interstellar turbulent pressure and 
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Fig. 1. — Surface area and HI masses: ADBS galaxies 
IIRosenberg Schneider 200j|) are plotted as black squares. The 
solid line is the theoretical curve fEguation llGp from this work and 
the area between the dashed lines represents the 2a variance region. 
Note the simple relation in which the HI mass scales almost lin- 
early with the surface area, implying a regulated HI surface density 
across galaxies. This hitherto unexplained feature is reproduced by 
our simple model. 

confirmed a very weak dependence of Ug on the star 
formation rate. They find simulated HI emission lines 
widths of 10 — 18 km s~^ for models with SN rates that 
range from 1 to 512 times the Galactic SN rate. Hence, 
the characteristic values for pd and ag when plugged into 
Equation[5l set the characteristic scale for the gas surface 
density to a ballpark figur e of ^ lO^Mtn kpc~^, which is 
close the observed value (jRosenberg fc Schneider! [200 3| ) 
for the ADBS galaxies. 

5. REGULATED SURFACE GAS DENSITY 

Integrating the st ar formatio n rate per unit volume 
over an entire galaxv lSilS (|1997[ ) provides the total SFR 
as 



A'h ^ 2.6 X i;. 



5/2 

rot, 200 



Q 



-3/2 



P 

(15 



-0.87 



Meyr- 



(9) 



where Wrot, 200 is the maximum rotational velocity of the 
galaxy, nor malized to 200 km s~^. The exact expression 
(jSilkl Il99"7[ ) depends on well understood quantities like 
the total mass of stars formed to yield one core collapse 
supernova, the mechanical energy output of an average 
supernova, and depends only weakly on the metallicity. 
As the disk becomes gravitationally unstable if Q de- 
creases, the SFR increases and the resulting supernovae 
increase the porosity P. This reduces the supply of cold 
gas a nd he nce suppresses the star formation rate. Hence 
ISilkl (I1997D points out that this expression provides an 
explicit demonstration of disk self-regulation and self- 
regulated star formation ensures P ^ 0.5 and Q ~ 1. 
We use these values in the rest of this work. 

Using the definitio n of the dynamica l mass {Md = 
{AV)'^Rgo%/G)) from lKulkarni fc HeilesI (IT988I) . we can 
express tVot, 200 as 



Vrot,200 ^ 1-04 X 



Md, 
R 



10 



(10) 



where Md,ii is the dynamical mass normalized to 
IO^^A/q and i?io is the 90% optical containment radius 
(i?go%) normalized to 10 kpc. The dynamical mass has 
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been found to be correlated with the cu be of the op- 
tical radius (iGarcia-Appadoo et al.l [20091 ). The implied 
roughly constant global dynamical density is of the order 
of pd ~ l(f M^kpc'^ . Exploiting this observed relation 
we have 

A/d,n ^ 0.42p7i??o- (11) 

where py ~ pd/ lO'^ MQkpc~^ . Substituting for Urot, 200 
and then for M^^ii, we get 

Af* i2 0.95 X 4i?fo Mq yr-^ (12) 

One may interpret the iGavazzi et al.l (|1996f ) relation be- 
tween total luminosity (a proxy for the stellar mass) 
and dynamical mass (proportional to because of the 
shared pd) as Af, cx R^. This implies that the doubling 
time for the stellar mass, scales as 

r=^cxi?i (13) 
Af, 

Objects with larger doubling times have older stellar pop- 
ulations on an average. This may explain why bigger 
galaxies are systematically redder. 

This star formation rate allows us to estimate the mean 
SFR per unit area as 

^SFR ^ 3.0 X 10-3 X p^-25i??o' Mq yr-1 kpc-^. (14) 

However ^sfr is related to the s urface dens i ty of gas 
Sgas by the Kennicutt-Schmidt Law lKennicuttI (|1998[ ) as 



Es™ =(2.5 ±0.7) X 10-4 

^ n n. .t; 



(15) 



1.4±0.15 



Mq yr-1 kpc- 



. 1 Mq pc"2 , 

If the molecular gas fraction is <C 1, which is true for 
most disk galaxies, the mean surface density of HI 'Shi 
is comparable to the total gas surface density Sgas- This 
assumption will make a proportional error comparable to 
the molecular gas fraction. Hence assuming E// / ~ Sgas , 
eliminating the 'Ssfr and integrating over the surface 
area, we obtain the relation between the gas mass and 
the surface area as 



log 



^ (6.91 ±0.08) + 1.18 X ;o.g 



lAf, 



A 



lkp<y 



- 0.89 X log 



Pd 



WMQkpc 



-3 



(16) 



where A = ttR^ is the cross sectional surface area pre- 
sented by the galaxy. Given that dynamical mass den- 
sity {pd) is not seen to vary much across galaxies, the 
almost linear relation between the total HI mass and 
the surface area implies very similar HI surface densities 
across a rang e of galaxies. The normalizat ion matches 
the observed ([Rosenberg fc Schneider] [20031 ) relation for 
P7 = 0.26 which lies within the range of its observed 
([Garcia-Appadoo et al.ll2009D values. The HI masses, of 
ADBS galaxies ([Rosenberg fc Schneidedl2003D . spanning 
3 orders of magnitude, vary only by ~ 0.14 dex (Icr) from 
the theoretical curve. However, this scatter is larger than 
the scatter propa gated from the Kennicutt-Schmidt Law 
([Kennicuttlll998| ). The excess scatter may be attributed 
to the scatter in core pr operties of halo s such as P7. It 
has been pointed out bv lLoeb fc Weiner[ ([2010f ) that nu- 
merical simulations are required to determine the scat- 
ter in dark matter core densities as a function of mass 



and rcdshift. Thi s would be important for dark matter 
dominated halos. IGavazzi et al.l ([1996D show a correla- 
tion between total luminosity and the dynamical mass. 
If most of the dynamical mass is provided by the stellar 
content, it would be important to study the scatter in 
this relationship. 

6. DISCUSSIONS 

Our result shows that the present understanding of 
mechanical feedback from supernova e, leading to self 
regulated star formation ([SilU |1997[) , can account for 
the regulation of HI surface density across galaxies 
to a characteristic value of around ~ l{f Mq kpc"^ 
([Rosenberg fc Schneideill2003| ). The predicted HI surface 
density depends on the mean density of the dynamical 
mass, which is likely to be provided by a combination 
of dark matter a nd stellar content. As to w hy galax- 
ies are observed ([Garcia-Appadoo et al.l [2009( 1 to share 
similar dynamical mass densities is still an open case re- 
quiring further investigation. Cold dark matter particles 
inter acting through a Yukawa potential (iLoeb fc Weineii 
|2010[ ) could provide a natural explanation for a charac- 
teristic den sity in dark matter dominated halos. The 
correlation ([Gavazzi et al.l[l"996[ ) between luminosity and 
dynamical mass may be important in halos dominated 
by the stellar content. 

The model of self regulated star formation ([SilU 1 19971 ) 
has been shown in this work to provide a scaling between 
doubling time and radius. This could explain why larger 
galaxies are systematically redder. This relation should 
be used in conjunction with population synt hesis mod- 
els such as Starburst99 ([Leitherer et al.iri999f ) to predict 
colors as a function of galaxy size. This could provide an 
interesting test of this model in future. 

Even if mergers drive a galaxy away from the funda- 
mental line, self regulation of the porosity of the ISM and 
Toomre parameter of the star forming disk can bring the 
HI surface density back to the value which is predicted 
in our simple model and observed in a wide range of 
galaxies. In a framework for hierarchical galaxy merg- 
ers this can happen multiple times in a galaxy's history, 
until it eventually runs out of neutral hydrogen. A de- 
tailed understanding of the proposed mechanism would 
require self consistent galaxy simulations taking into ac- 
count the supply of hot gas into the ISM from individual 
supernovae. 

Our model relates the total neutral hydrogen mass of 
the galaxy with its projected surface area. However in 
practice, the quantities observed with radio telescopes 
are the redshift, integrated HI line fluxes and the solid an- 
gles on the sky. Fluxes and solid angles behave differently 
in different cosmological scenarios, as they scale with the 
luminosity distance and angular diameter distance cIa 
respectively. This could facilitate their use in testing the 
lEtheringtonl ([1933( 1 relation, = (1 + z^dA, when fu- 
ture telescopes such as the Square Kilometer Array start 
to detect red-shifted HI from very distant galaxies. 
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